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(§4) PROCEDE DE FILTRAGE NUMERIQUE ADAPTATIF DANS LE DOMAINE FREQUENTIEL. 

fe7; L'invention concerne un procede de filtrage adaptatif 
aesignaux d'entree (x(n)) temporels dans le domaine fre- 
quentiel. Le filtrage est realise sur des blocs (xjk)) de m si- 
gnaux successrfs ((x(n), x(n+m-1))), par multiplication 
dans le domaine frequentiel de la transformee de Fourier 
(X (k)) des blocs de signaux d'entree par des blocs 
(W^(k)) de coefficients de ponde ration, ces blocs de coeffi- 
cients etant calcules recursivement a partir de blocs de 
coefficients precedents et de premier et deuxieme termes 
correctrfs (C^k), C 2 (k)) obtenus d'une part a partir de si- 
gnaux de sortie (y m (i<)) correspondents aux signaux d'en- 
tree filtres et d'autre part a partir de la transformee de 
Fourier d'un modele de I'entree (b^). 
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PROCfiDfi DE FILTRAGE NUMERIQUE ADAPTATIF 
DANS LE DOMAINE FREQUENTI EL 

L f invention concerne un procede de filtrage 
numerique adaptatif, dans le domaine f requentiel , de 
signaux numeriques. De tels procedes sont mis en oeuvre 
en traitement du signal, pour adapter les dispositifs 
5 de reception aux evolutions d'un canal de transmission 
evoluant temporellement . Ainsi, par exemple, en 
transmission hertzienne, les caracteristiques de 
propagation des signaux evoluent en fonction des 
conditions meteorologiques : elles sont differentes dans 
10 le temps et dans l'espace. On utilise done 
generalement, dans les recepteurs, des filtres 
adaptatif s, e'est a dire des filtres a coefficients de 
ponderation variables dans le temps- Les variations 
dans le temps des coefficients de ponderation sont 
15 definies selon un critere d • opt imisation et ces 
coefficients sont produits par des dispositifs 
(generalement des processeurs de traitement de signal 
ou DSP) mettant en oeuvre un algorithme d 1 adaptation . 

Des coefficients de ponderation, calcules par des 
20 algorithmes d' adaptation generalement employes, 
peuvent, dans le domaine temporel, se mettre sous la 
forme iterative suivante (les indices representant les 
dimensions de vecteur ou de matrice) : 

w m (n) = w m (n-l) - MH m - 1 {J(n) } .G{J(n) } (1) 
25 avec : 

n indice temporel, 

m 1 1 ordre du filtre, 

fi le pas de convergence, 

J(n) = E{|d(n) - Y(n)| 2 } (fonction de cout 
30 minimisee) , 

y(n) les signaux de sortie de filtre correspondants 



2738692 

2 

a des signaux d'entree *(n), 

d(n) les decisions produites par un organe de 
decision a partir des signaux de sortie, 

H m ^ 1 {J(n)} une estimation de 1 1 inverse de I'Hessian 
5 H de J(n) , et est une matrice de dimension m x m. 

et G{J(n)} une estimation du gradient de J (n) et 
est un vecteur de dimension m. 

Cette formulation est tres generale et est valable 
pour une tres grande classe d 1 algor ithmes echantillon 
10 par echantillon dans le domaine temporel. Les 
algor ithmes correspondant a 1» equation (1) sont dits de 
type quasi-Newton. 

Un choix connu de G{J(n)} est : 

G{J(n)} = - 2 x m *(n) .e(n) , avec 
15 x m( n * = t x ( n )#— -/ x ( n - m + 1) ] T un vecteur forme 

des m dernier signaux successifs recus au temps n, * 
etant le symbole de la conjugaison complexe, T etant 
une transposition ligne-colonne , et 

e(n) = d(n) - x m T(n) - w m (n " x) un si< ?nal d'erreur. 
20 C'est le choix de H m ^ 1 {J(n)} qui caracterise 

principalement les differents algorithmes de quasi- 
Newton . 

Si cette matrice est la matrice identite, on 
obtient l'algorithme des moindres carres moyens (ou 

25 LMS, de l'anglais least mean squares). Cet algorithme, 
tres utilise, presente une complexite faible, au 
detriment d'une convergence lente egalement. 

Si cette matrice est l f inverse de la matrice 
d • autocorrelation des signaux d'entree, alors on 

30 obtient l*algorithme des moindres carres recursif (ou 
RLS, de l'anglais recursive least squares). Cet 
algorithme presente une convergence meilleure que 
1' algorithme LMS (les signaux d'entree etant pris en 
compte) , au prwix d'une complexite beaucoup plus 
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importante. Un probleme de mise en oeuvre de cet 
algorithme d 1 adaptation est le temps necessaire a la 
realisation du calculs de H^" 1 et G pour chaque bloc de 
signaux d' entree. 

D'autres algorithmes de filtrage ont ete proposes, 
derives du type RLS, bases sur une extrapolation min- 
max d'une matrice d • autocorrelation d'un ordre p + 1, 
inferieur a m, jusqu'a 1 ■ ordre m du filtre (G.V. 
Moustakides, S. Theodoridis, "Fast Quasi Newton 
Transversal Filters - A New Class of Adaptive 
Estimation Algorithms", IEEE Trans. on Signal 
Processing, vol. 39, pp 2184-2193, Oct. 1991). 

On utilise alors une estimee de 1 1 inverse de la 
matrice Hessienne donnee par la formule suivante : 



R m" 1 (") = 



R p -1 (n) 0 



+ r. 



a_-!(n - j) 



b m b m H ( 2 ) 
mm x ' 



1, avec p < m, 
de dimension 



J D "P 

0 0 

avec 

j indice variant de 0 a m - p - 
Rp" 1 ^) une matrice carree 
completee de m-p zeros, 

b m " [OjTr-bpHcn - j ) : 1 : 0 m _ p _ -5.^] H , 
b m H la transposee hermit ienne de b m , 

b p (t) et a p (t) les predicteur recursif et puissance 
d'erreur de prediction d' ordre p, 

l'operateur : representant une concatenation par 
juxtaposition, 0j T etant un vecteur nul de dimension j 
transpose, 0 m- p -j-i T ^ tant un vecteur transpose d'un 
vecteur nul de dimension m-p-j-1. 

lis permettent d'atteindre une convergence 
equivalente a un algorithme de type RLS tout en ayant 
une complexity, proche de celle d'un algorithme de type 
LMS dans le domaine temporel, necessitant 2m + 6p 
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multiplications par unite de temps. De tels procedes 
sont neanmoins de complexite elevee pour des filtres 
d ( ordre eleve, du fait de leur complexite 
proportionnelle a m (ce qui peut rendre prohibitive 
5 leur utilisation pour des applications telles que 
l'annulation d'echos en audioconf erence , ou m peut etre 
de l'ordre de plusieurs milliers) . 

Depuis quelques annees, dans le domaine du 
traitement du signal, le filtrage adaptatif dans le 

10 domaine frequentiel a fait I'objet de nombreux 
developpements . Le filtrage adaptatif dans le domaine 
frequentiel presente l'avantage, par rapport au 
filtrage adaptatif dans le domaine temporel, de reduire 
considerablement la complexite de la mise en oeuvre des 

15 algorithmes d 'adaptation. Ceci provient du fait que la 
convolution temporelle est remplacee dans le domaine 
frequentiel par une multiplication et deux transformees 
de Fourier. En consequence, un calcul dans le domaine 
frequentiel permet de reduire a la fois les temps de 

20 production des coefficients de ponderation et 
1 • encombrement des circuits calculant ces coefficients. 

Les algorithmes d' adaptation dans le domaine 
frequentiel mettent en oeuvre un filtrage adaptatif par 
blocs, avec une dimension de bloc generalement egale a 

25 Pordre du filtre, c'est a dire au nombre de 
coefficients de ponderation du filtre. La plupart des 
algorithmes existants sont du type algorithme du 
gradient. On connait deja des exemples de mise en 
oeuvre d 1 algorithmes LMS dans le domaine frequentiel. 

30 Ainsi, E.R. Ferrara, "Frequency Domain Adaptive 
Filtering", Adaptive Filters, C. F. N. Cowan and P. M. 
Grant, Eds, Englewood Cliffs, Prentice Hall, 1985, 
Chapitre 6, page 145-179 propose des filtres adaptatif s 
dans le domaine frequentiel. Plus particulierement , la 
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figure 6.2 et la description correspondante presentent 
un filtre adaptatif de type LMS dans le domaine 
frequentiel. En terme de complexite de realisation, le 
calcul de m echantillons de sortie d'un filtre a m 
5 coefficients avec 1 1 algorithme LMS dans le domaine 
temporel requiert 2m 2 multiplications reelles. Avec 
1 • algorithme decrit dans ce document le meme calcul 
necessite cinq transformees de Fourier rapides (encore 
notees FFT) de dimension 2m et deux multiplications 

lO complexes sur 2m points. Une FFT sur 2m points demande 
(m/2) log 2 m - m multiplications complexes. Le nombre de 
multiplications complexes par bloc donne est (5m/2) 
log 2 m. Une multiplication complexe equivalent a quatre 
multiplications reelles, le rapport de complexite entre 

15 LMS frequentiel et LMS temporel en terme de 
multiplications reelles est de (5 log 2 m + 4)/m. Pour m 
= 512 par exemple, le rapport sera de 0.096 (ce qui 
illustre 1 • interet de la mise en oeuvre dans le domaine 
frequentiel en terme de complexite de calcul) . 

20 Un but de 1" invention est de proposer un procede de 

filtrage adaptatif, dans le domaine frequentiel, ayant 
une convergence plus elevee que le filtrage LMS par 
blocs dans le domaine frequentiel. Pour ce faire, 
l 1 invention propose de modifier le procede de calcul 

25 des blocs de coefficients de ponderation en prenant une 
estimation du Hessian differente de l'identite (c'est a 
dire a produire a partir des signaux d 1 entree un modele 
de I'entree, ce qui permet d'ameliorer la convergence), 
et qui ne necessite pas le traitement de recursivite 

30 inherent a la formulation de type RLS . 

Plus precisement, 1' invention propose de transposer 
dans le domaine frequentiel une version simplifiee du 
procede illustre ci-dessus par l'equation (2). 

L 1 invention est basee sur la supposition que la 
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matrice d • autocorrelation d 1 entree ne change 
pratiquement pas pendant un intervalle de temps 
equivalent a au moins m echanti lions, c^st a dire 
qu'une estimee de la matrice d • autocorrelation d'entree 
qui est exacte a un temps n sera exacte pour au moins m 
a 2m intervalles de temps plus tard (un intervalle de 
temps correspondant a la distance temporelle entre deux 
signaux d'entree successifs) . De nombreuses 
applications repondent a cette hypothese pour des 
intervalles de temps beaucoup plus importants (par 
exemple l'annulation d'echos induit par une propagation 
en trajets multiples dans la transmission hertzienne) . 

Si on considere que ^(n) est constante pour au 
moins m ou 2m intervalles de temps, alors R p (n) le sera 
aussi • 

On peut done simplifier ^equation (2) de telle 
sorte que pour j variant de 0 a (m-1) - p on a : 

a p (n - j) = a p (n) et b p (n - j) = b p (n). 

En considerant que G{J(n)} = - 2 x m *(n) e(n), et en 
ecrivant liquation (1) m fois, avec n = km,..., km + m 
-1, on obtient une nouvelle equation recursive entre 
une estimation des coefficients de ponderation et 
1* estimation des coefficients de ponderation consideree 
m intervalles de temps plus tot : 

w m (k+l) = w m (k) + 2m £j i^Vm^.jW ( 3 >' avec 

j variant de km a km + m - 1, 

j (k) les elements du bloc d 1 erreur defini par : 

*m'( k > = d m< k > " x m T ( k ) w m( k ) < 4 >' et 

X m (k) = [x m (km) x m (km+m-l) ] vecteur d'entree 

k etant le rang du bloc. 

En supposant que la matrice d • autocorrelation 
d'entree est invariante pour m intervalles de temps 
successifs, on en deduit que : 

w m (k + 1) - w m (k) + 2m R rn - 1 X m *(k)e ni (k) (5), 
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la matrice R m ^ 1 (k) etant celle correspondant au 
dernier instant (n = km + m - 1) du k-ieme bloc. 
On peut encore ecrire : 



R^^k) = P m (k) + B(k).B H (k) avec 
Rp^ 1 (k) o 

0 0 



P m (k) = 



10 



15 



20 



25 



30 



et B m (k) est une matrice m * p definie par 



B m (k) = 



^p+l 



0 r 



'm-p-1 w m-p-2 



'm-p-l 
b p+l 



avec b p+1 = a p -l/2(k) [-b p H(k) 1]H. 
On peut ainsi mettre 1' equation (3) sous la forme : 
w m (k + 1) = w m (k) + 2 M c m l(k) + 2m c m 2 (k) (6) avec 

= P m( k )-V< k >- e m< k > < 7 >< et 
Cin 2 (k) = B (k) .B H (k) -X m *(k) -e m (k) (8). 

On notera que seuls les p premiers elements de 

c m 1 ( k ) sont non nuls. Pour le calcul du vecteur c m 1 (k) 

une inversion matricielle et une multiplication 

matricielle sont necessaires pour chaque bloc. En 

choisissant p << m on pourra choisir : 



"TTl 

matrice 
^(k). 



[ (X *(k) 



e m (k)) 



H 



Hi H 



avec 



p.m comprenant les 



• 0 "1 
p premieres 



X p *(k) 
rangees 



une 
de 

de 
au 



Autrement dit, le premier terme correct if 
1 'equation (6) est traite de maniere analogue 
procede LMS. Si on fait abstraction du terme correctif 
c m 2 (k) dans ^equation (6), on retrouve une equation de 
type LMS (a la seule difference que la contrainte de 
gradient dans le cas de c m 1 (k) est differente du fait 
que seuls p elements de ce vecteur sont non nuls) . 
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Les convolutions lineaires des equations (4) pour 
produire le bloc d'erreur et (7) pour le premier bloc 
correctif peuvent etre aisement implantees dans le 
domaine frequentiel en utilisant la technique de 
recouvrement partiel / dite overlap-save, bien connue de 
l'homme du metier. L'idee de base est de traiter les 
convolutions lineaires par le biais de convolutions 
circulaires de taille double, ces convolutions etant 
mises en oeuvre en utilisant des transformees de 
Fourier discrete rapide. Ainsi, on obtient : 
c 2n X M = FQ 1 F-lx 2m H(k) .E 2m (k) (9) avec 
E 2n| (k) le vecteur bloc d'erreur dans le domaine 
frequentiel, 

F et F"" 1 l'operateur de transformee directe de 
Fourier et son inverse, et 

q! une matrice de dimension 2m x 2m telle que : 




o 



0 0 

avec I la matrice identite de dimension m 
m 

complete, pour Q 1 , par des zeros. 

Si on regarde le second terme correctif de 
^equation (6), on peut verifier d'apres l'equation (8) 
que ce terme a une structure particuliere et peut etre 
ecrit sous forme d'une succession de trois convolutions 
lineaires . 

En terme frequentiel on pourra se ramener a la 
formulation suivante : 

C 2m 2 (k) = FQ2F-lB 2m *(k)FQ3F-lB 2m (k)FQ2 F - 
lx 2m H (k)E 2ln (k) (10) avec 

B 2m (k) definie par : 

B 2jn (k) = diag{FFT[b* p+11 :O in _ p _ 1 T :b* p+lp+1 
b* p+12 ) T } avec b* p+li le i-ieme element du vecteur 

b *p+l' 

q2 e t q3 une matrice de dimension 2m x 2m telle 
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que: 



Q 2 = 



-m 



et 



)3 = 



avec I m _p la matrice identite de dimension m - p. 

On remarque que 1' equation (10) implique cinq 
transformees de Fourier supplementaires par rapport a 
la transposition de l'algorithme LMS dans le domaine 
10 frequentiel, ces transformees etant utilisees pour 
mettre en oeuvre dans le domaine frequentiel des 
contraintes temporelles . 

Si on considere que p << m, alors on peut 
simplifier 1» expression de C 2m 2 (k) de sorte que : 
15 C 2ln 2 (k) = B* 2Tn (k).B 2m (k).X 2m H (k) .E 2Tn(k) 

Cette simplification est raisonnable parce que, en 
supposant que p est tres inferieur a m, alors les 
effets de distorsion sont negligeable et les 
contraintes matricielles peuvent etre supprimees. 
20 Si on considere, dans le domaine frequentiel, 

l'equation recursive d'adaptation des coefficients de 
ponderation suivante : 

W 0 «(k+1) = W 2ln (k) + 2 M C 2l /(k) + 2 M C 2TO 2(k) (12) avec 
et C 2m 2 (k) tels que definis par les relations 
25 (9) et (11), alors on obtient un procede de filtrage 
dans le domaine frequentiel ayant une complexity proche 
de celle de la transposition f requentielle de 
l'algorithme LMS, et presentant une convergence 
analogue a celle de 1 1 implantation dans le domaine 
30 temporel de l^lgorithme RLS defini dans la relation 
(5) . 

Ainsi, !• invention concerne un procede de filtrage 
adaptatif de signaux d 1 entree temporels, caracterise en 
ce que le filtrage est realise sur des blocs de m 



"2m 
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signaux successifs, par multiplication dans le domaine 
frequentiel de la transformee de Fourier des blocs de 
signaux d' entree par des blocs de coefficients de 
ponderation, ces blocs de coefficients etant calcules 
recursivement a partir de blocs de coefficients 
precedents et de premier et deuxieme termes correct if s 
obtenus a partir de signaux de sortie correspondant aux 
signaux d 1 entree filtres et a partir de la transformee 
de Fourier d f un modele de 1' entree. 

D'autres particular ites et avantages apparaitront a 
la lecture de la description qui suit, a lire 
conjointement au dessin annexe dans lequel la figure 1 
represente schema tiquement un procede de filtrage 
adaptatif dans le domaine frequentiel, selon 
1 1 invention. 

La figure 1 illustre schematiquement un exemple de 
mise en oeuvre du procede de filtrage adaptatif selon 
1 ' invention. 

L'objet du procede est de ponderer des signaux 
d'entree x(n) successifs par des coefficients w i (k) 
variables et de produire des signaux de sortie y(n) 
correspondants tels qu'on minimise des signaux d'erreur 
e(n) = d(n) - y(n) r avec d(n) des decisions produites 
par un organe de decision a partir des signaux de 
sortie y (n) . 

Une application d 1 un tel procede est par exemple 
l'annulation d'echos en audioconf erence . Une autre 
application possible est 1 • identification adaptative 
d'une maniere generale, et en particulier 

1 • identification de canaux perturbes par des echos lors 
de transmissions ou de diffusions hertziennes. Dans ce 
contexte, l'ordre des predicteurs ainsi que les 
conditions de stationnar ite rentrent exactement dans 
les hypotheses de 1' invention. 
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Dans 1' invention, les signaux d' entree x(n) sont 
traites par blocs de m signaux successifs (x(n),..., 
x(n+m-l)), avec m entier. Generalement m est une 
puissance de 2 (par exemple m = 512 ou 1024). 
5 Class iquement, a un bloc de signaux d * entree 

(x(n),..., x(n+m-l)) va correspondre un bloc de m 
coefficients de ponderation (w n (k),..., w n+m _ 1 (k)). 

Schema t iquement, le procede selon 1* invention 
propose de realiser, dans le domaine frequentiel, 

10 l'adaptation des coefficients d'une part a partir d'un 
bloc d'erreurs produit a partir des signaux de sortie 
(de roaniere analogue a la realisation dans le domaine 
frequentiel d f un procede de type LMS) et d' autre part a 
partir d'un predicteur de signaux d» entree. 

15 on va decrire maintenant plus precisement le 

procede illustre sur la figure. 

Les signaux d'entree x(n) sont regroupes, par une 
transformation serie-parallele 1, en blocs de signaux 
(x(km) , . . . , x(km+ra-l) ) . 

20 Chaque k-ieme bloc x^Ck) est regroupe par une 

concatenation 2, avec le bloc precedent x^k-l), puis 
le bloc de 2m signaux obtenus est transpose dans le 
domaine frequentiel, par une transformee de Fourier 
rapide 3 (notee FFT) . On obtient alors un bloc 

25 frequentiel de 2m signaux d 1 entree. 

Par une multiplication 4, on produit 2m signaux 
representatif s du resultat de 1' operation suivante : 

Y 2m< k > - X 2m( k > W 2m( k ) avec 

X 2m (k) = diag{FFT[x(km-m) . . . .x(km+m-l) ] T } et 

30 rk) un vecteur, de dimension 2m, de coefficients 

^ 2 m x 1 

de ponderation. 

Apres une transformation de Fourier rapide inverse 
5 de Y 2Tn (k) et une deconcatenat ion 6 on produit un bloc 
y (k) de m signaux correspondant a la seconde moitie 
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des signaux obtenus apres la transformation inverse, 

dans le domaine temporel. 

Autrement dit, le bloc y m (k) est representatif de 

la relation suivante : 

5 y m (k) = Q°F- 1 Y 2m (k) avec 

qO - r 0 I m ] une matrice de dimension m.2m, avec 
L in * m j 

0 m une matrice carree nulle de dimension m, I ro la 
matrice identite de dimension m, et F"" 1 1 1 inverse de 
l'operateur de transformee directe de Fourier de 

10 dimension 2m. 

Par une transformation parallele-serie 7, on 
obtient les signaux de sortie y(n). 

En ce qui concerne 1 'adaptation des coefficients de 
ponderation, on realise 1' operation suivante : 

15 W 2m (k+1) = W 2Tn ()c) + 2MC 2m l(k) + 2MC 2m 2 (k) (12 ), 

avec 

C 2m l(k) = FQlF-lx 2m H(k).E 2ln (k) (9) et 

c 2m 2 < k > = B 2m*( k >- B 2m( k >- x 2m H < k >- E 2m< k > 
avec 

20 m un pas de convergence, 

E 2m (k) un vecteur bloc d ■ erreur dans le domaine 
frequentiel, de dimension 2m, 

X 2m H (k) une matrice diagonale de dimension 2m. 2m 
obtenue par une transposition hermitienne 8 de X 2m (k) , 
25 F l'operateur de transformee directe de Fourier de 

dimension 2m, 



Q 1 une matrice de dimension 2m * 2m telle que : 



30 




p entier inferieur a m, 

Ip matrice identite de dimension p, 

B<>_(k) matrice de dimension 2m definie par : 



35 
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B 2m< k > = dia g< FFT C b *p+l / l :0 TT»-p-l T:b *p+l / p+l 

b* , , ^1 T > avec b* ... ; le i-ieme element d • un vecteur 
p+l,2 J J p+1,1 

b*p +1 de dimension p+1. 

Pour produire ^^(k), on P roc ^ de de la maniere 
5 suivante : 

- a partir du bloc y m (k) de m signaux de sortie, 
definis dans le domaine temporel, et d'un bloc d m (k) de 
m decisions correspondantes , on produit par une 
soustraction 9 un bloc de m signaux d'erreur ^(k) = 

10 d m( k ) " Y m ( k ) • Les decisions sont produites par un 
circuit de decision non lineaire, 

- apres une concatenation 10 de e m (k) avec un bloc 
de m signaux nuls, on produit le bloc E 2m (k) de 2m 
signaux d'erreurs dans le domaine frequentiel par une 

15 transformee de Fourier rapide 11 (on a done E 2m (k) = 
FQ 0T e m (k)), 

a partir de E 2m* k ) on produit, par une 
multiplication 12 avec le bloc X 2m H (k) de 2m signaux, 
un bloc de 2m signaux representatif s de I'operation 
20 X 2ln H (k) . E2m(k)/ 

ce bloc de 2m signaux obtenu est ensuite 
transpose dans le domaine temporel par transformee de 
Fourier inverse rapide 13, les m derniers signaux 
correspondants sont remplaces par des zeros lors d'une 
25 concatenation 14, et les 2m signaux ainsi produits sont 
transposes dans le domaine frequentiel par une 
transformee de Fourier rapide 15. 

Autrement dit, le terme correctif C 2m 1 (k) est 
obtenu par un procede semblable a celui utilise dans la 
30 mise en oeuvre d 1 un algorithme de type LMS dans le 
domaine frequentiel, en remplacant I 2m par Q 1 . 

Le terme correctif C 2m 2 (k) correspond a une 
modelisation des signaux d r entree. 

L 1 invention est basee sur la supposition que la 
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matrice d 1 autocorrelation d' entree ne change 
pratiquement pas pendant un intervalle de temps 
equivalent a au moins m echantillons, c'est a dire 
qu'une estimee de la matrice d ' autocorrelation d 1 entree 
5 qui est exacte a un temps n sera exacte pour au moins m 
a 2m intervalles de temps plus tard (un intervalle de 
temps correspondant a la duree entre deux signaux 
d f entree successifs) r voire plus. De nombreuses 
applications repondent a cette hypothese pour nombre 
10 l.m d 1 intervalles de temps beaucoup plus important (par 
exemple l'annulation d'echos induits par une 
propagation en trajets multiples dans la transmission 
hertzienne) . 

En consequence, on considere des regroupements dans 

15 un concatenateur 16 de 1 blocs successifs de signaux 
d 1 entree. Pour chaque regroupement, on calcule dans le 
domaine temporel un bloc b p+1 de p+1 signaux (modele de 
1* entree) tels que 

b p+i = a p - 1 / 2 (k) [-b p H(k) 1] H avec 

20 bp et a p des predicteur d* entree et puissance 

d'erreur de prediction d'ordre p, avec p < m. 

Pour ce faire, on utilise dans un correlateur 17 la 
methode connue dite d • autocorrelation (decrite par 
exemple dans "Digital Spectral Analysis with 

25 Applications", par S.L. Marple, Prentice Hall, New 
Jersey, 1987) . 

A partir du bloc b p+i on produit par une 
concatenation 18 dans le domaine temporel le bloc b' m 
de m signaux (b* p+1 ± 0 m _ p _ 1 T b* +1 p+1 . . . b* p+1 2 ) - 

30 On passe ensuite dans le domaine frequentiel par 

une transformee de Fourier rapide 19 pour obtenir un 
bloc de 2m echantillons representatif de B 2m (k) tel 
que: 

B 2m Ck)-diag<FFT(b- p+1#1 sO |n _ p _ 1 T:b- p+liP+1 -..b* p+1#2 lT>. 
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Pour obtenir C 2in 2 (k) , on precede ensuite a la 
multiplication 20 de B 2in (k) par sa conjuguee B 2m *(k) 
(obtenue par une transformation 21), et a la 
multiplication 22 par X 2m H (k) . E 2m (k) , ce terme ayant 
5 ete calcule par ailleurs pour produire C 2ln 1 (k). 

A partir des deux blocs C 2in 1 (k) et C 2m 2 (k) de 2m 
signaux correctifs, et apres multiplications 22 et 23 
de ces blocs par 2p, on peut calculer pour chaque bloc 
de signaux d' entree le bloc de coefficients de 
10 ponderation correspondant par une simple addition 24 
des signaux du bloc de coefficients de ponderation 
precedent fourni par un circuit de retard 25 avec les 
signaux des blocs 2MC 2m 1 (k) et 2nC 2m 2 (k) . 
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REV END I CATIONS 

1 - Procede de filtrage adaptatif de signaux 
d'entree (x(n)) temporels, caracterise en ce que le 
filtrage est realise sur des blocs (x m (k)) de m signaux 
successifs ( (x(n) , . . . ,x(n+m-l) ) ) , par multiplication 
5 (4) dans le domaine frequentiel de la transformee de 
Fourier (X 2in (fc)) des blocs de signaux d • entree par des 
blocs ( w 2m^ k ^ de coefficients de ponderation, ces 
blocs de coefficients etant calcules recursivement (9- 
15) a partir de blocs de coefficients precedents et de 

10 premier et deuxieme termes correctifs ( c 2m 1 ^ k) ' 
C 2m 2 (k)) obtenus d'une part (13-15) a partir de signaux 
de sortie (y m (k)) correspondants aux signaux d'entree 
filtres et d' autre part (16-20) a partir de la 
transformee de Fourier d 1 un modele de 1 • entree (b p+1 ) • 

15 2 - Procede selon la revendication 1, caracterise 

en ce que le modele de 1' entree (b p+1 ) est calcule par 
autocorrelation (17) et comprend un nombre p de 
signaux, avec p infer ieur a m. 

3 - Procede selon la revendication 2, caracterise 
20 en ce qu'un modele de 1 1 entree est calcule (16) pour 

l.m blocs de signaux d 'entree successifs. 

4 - Procede selon l'une des revendicat ions 1 a 3, 
caracterisee en ce que le calcul du second terme 
correctif ( c 2m 2 ^ k ^ > comprend la transformation par 

25 transformee de Fourier (19) d 1 un bloc (^* m ) de m 
signaux produits dans le domaine temporel a partir du 
modele de 1 1 entree (b p+1 ). 

5 - Procede selon l'une des revendications 1 a 4, 
caracterise en ce que le premier terme correctif 

30 ( c 2m 1 ^ k ^ est obtenu P ar la mise en oeuvre (10-15) dans 
le domaine frequentiel de 1 ' algor ithme des moindres 
carres moyens . 
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